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HARD  X-RAY  SPECTRUM  FROM  A  ONE-FLUID  LASER  PLASMA  MODEL 
I.  Introduction 

The  space-time  integrated  experimental  x-ray  spectra  obtained  from 
laser-produced  plasmas  typically  show  a  two-temperature  structure.  The 
energetic  part  of  the  spectrum  -  which  yields  the  hot  temperature  T^  -  has 

been  used  as  supporting  evidence  for  suprathermal  electrons^,  flux  inhi- 
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bition  mechanism  and  self -generated  magnetic  fields.  The  first  inter¬ 
pretation  is  however  the  most  common  one  and  relies  on  the  generation  of 
a  small  suprathermal  electron  population  directly  responsible  for  the 
energetic  X  rays.  This  population  is  assumed  to  coexist  everywhere  with 
the  bulk  of  the  electrons  which  are  much  cooler  and  which  give  rise  to  the 
cold  temperature  T£  in  the  x-ray  spectrum. 

Generally  in  presenting  hard  x-ray  data,  lines  are  drawn  through 
experimental  points  which  show  sharp  breaks  in  the  curve.  These  breaks 
are  used  to  justify  the  coexistence  of  two  electron  temperatures  every¬ 
where.  However,  we  should  note  that  a)  some  of  these  kinks  occur  even  at 
relatively  small  energy  (below  10  keV)  and  b)  often  smooth  curves  can  be 
drawn  through  the  error  bars  which  fit  the  data  as  well  as  kinked  curves. 

We  show  on  various  specific  cases  that  in  accordance  with  ref.  2 
and  3  for  example,  it  is  also  possible  to  interpret  the  data  in  terms  of 
a  single-temperature  electron  fluid.  One  of  the  differences  between  a 
one-temperature  and  a  multi-group  treatment  of  the  electrons  is  that  it 
is  easier  to  fit  an  experimental  curve  with  a  multi-group  treatment  since 
Maanaeript  wbmitUd  April  10, 1980 


several  parameters  (for  Instance  temperature,  number  of  hot  electrons  and 
flux  limit  for  each  energy  bin)  have  to  be  included  into  the  multi-group 
model.  In  a  fluid  model,  all  the  electrons  become  very  hot  for  a  short 
period  of  time,  in  a  very  localized  region.  The  basic  criterion  for  this 
model  to  hold  is  that  the  fluid  approach  has  to  be  valid.  In  that  case, 
the  mean  free  path  of  an  electron  or  its  Larmor  radius  must  be  smaller  than 
a  characteristic  length  of  the  plasma  (henceforth,  we  will  not  consider 
magnetic  fields).  For  a  Nd  laser-produced  plasma,  a  characteristic  length 
near  the  critical  density  is  of  order  ly.  At  other  places  in  the  plasma, 
it  is  at  least  10  times  as  long.  The  mean  free  path  of  a  10  keV  electron 
at  the  critical  density  is  about  0.2  cm.  Thus,  for  a  single  fluid  theory 
to  be  valid,  there  must  be  some  anomalous  process  which  greatly  reduces 
this  mean  free  path.  Our  model  invokes  return  current  driven  ion  acoustic 
turbulence  to  reduce  the  mean  free  path. 

In  this  work,  we  present  calculated  x-ray  spectra  from  a  one-temperature 
fluid  model.  We  have  calculated  these  spectra  for  different  energies  and 
we  compare  them  to  detailed  experimental  spectra  obtained  at  the  Naval 
Research  Laboratory.^  Several  variations  of  parameters  are  investigated 
and  discussed. 

II.  Model 

The  model  which  has  been  used  in  this  study  has  been  described  in 
detail  before. ^  It  is  characterized  by  three  different  absorption  mechanisms, 
stimulated  Brillouln  backscattering  and  heat  flux  limitation.  The  three 
absorption  processes  are  inverse  bremsstrahlung,  resonant  absorption  and 
absorption  due  to  ion  acoustic  turbulence.  This  latter  process  allows  the 
fluid  description  because  of  the  anomalous  collision  frequency.  In  ref.  6, 


it  has  been  shown  that  the  mean  free  path  is  given  roughly  by  V#/w  (^  ) 
where  (-p1)  is  a  measure  of  the  ion  acoustic  fluctuation  level.  If 

e  _i 

(e$/Tg)  ~  10  ,  the  electron  mean  free  path  for  a  10  keV  plasma  at 

21  -3 

n  *  10  cm  is  roughly  2p,  which  is  comparable  to  the  shortest  charac¬ 
teristic  length  in  a  Nd-laser  produced  plasma.  Heat  flux  limitation 
follows  from  the  same  physical  mechanism  and  is  included  in  a  self- 
consistent  fashion  into  the  code.  Stimulated  Brillouin  backscatter  is 
also  included  into  the  code  but  does  not  play  a  major  role  in  the  short 
pulses  considered  in  this  work. 

One  difference  between  this  work  and  ref.  6  is  the  treatment  of  the 
dense  region  behind  the  absorption  region.  In  the  present  studies, 
appreciable  amounts  of  radiation  come  from  this  region.  Therefore  more 
care  was  taken  in  modeling  this  part  of  the  plasma.  Specifically,  many 
more  grid  cells  were  included  in  the  high-density  region.  Also  the 
boundary  condition  at  the  high  density  side  was  taken  to  mock  up  thermal 
contact  with  a  slab  at  Tfi  ■  100  eV  and  also  having  -|^  ■  0,  -|^  ■  0,  to 
model  free  flow  of  matter  and  momentum  into  or  out  of  the  slab.  Since  a 
one-dimensional  model  might  be  inappropriate,  we  have  crudely  investigated 
2-0  effects  by  assuming  that  lateral  thermal  conduction  can  be  modeled  by 
a  loss  term  -K/Lc2  Te  in  the  temperature  equation  where  L£  is  a  thermal 
conduction  length  and  K  is  the  thermal  conduction  coefficient.  The  effects 
of  this  length  on  the  results  are  shown  in  detail  for  one  case.  The  free- 
free  radiation  was  broken  into  17  frequency  groups  going  from  700  eV  to 
100  keV.  These  groups  were  not  divided  evenly  over  the  whole  frequency 
range  in  order  to  provide  better  resolution  at  the  low  energy  range  of  the 
spectrum.  The  x-ray  emission  is  displayed  in  erg/keV  and  corresponds  to 
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radiation  emitted  in  4ir  steradians.  The  target  is  typically  a  CH  planar 
target  with  a  lOu  initial  characteristic  gradient  length  in  front  of  it 
(below  critical  density).  The  laser  pulse  is  gaussian  with  SO  psec  FWHM 
and  a  focal  spot  radius  of  15p  (the  focal  spot  radius  is  needed  in  order 
to  yield  absolute  radiation  levels). 

III.  Results 

We  attempted  to  reproduce  a  series  of  experiments  performed  at  the 

Naval  Research  Laboratory^  in  a  systematic  fashion.  This  work  is  then 

different  from  any  other  published  work  involving  one-temperature  models 

because  it  deals  with  a  whole  range  of  experiments.  The  experimental  data'* 

was  taken  between  2  x  If/**  and  8  x  10^  W/cm2  with  a  hard  x-ray  tail 
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appearing  above  7  x  10  W/cm  .  The  numerical  results  have  been  obtained 

15  2 

for  the  range  1.5  to  8  x  10  W/cm  for  4  values  of  the  irradiance:  1.5  - 
2.5  -  5  and  8  x  1015  W/cm2. 

No  prepulse  was  present  in  the  experiment  and  so  a  lOy  initial 
density  gradient  was  assumed  in  front  of  the  target  as  in  the  case  of  the 
short  pulse  simulations  in  ref.  6.  Because  of  this  sharp  density  gradient, 
backseat ter ing  is  not  significant. 

Before  presenting  results  for  the  x-ray  spectra,  let  us  note  that  the 
calculated  total  absorption  varied  from  a  high  432  for  the  1.5  J  case  down 
to  322  for  the  8  J  case.  Resonant  absorption  respectively  accounted  for 
272  and  232  for  these  two  cases. 

Results  and  comparisons  for  the  cases  mentioned  appear  in  Fig.  1-4. 

The  dotted  line  for  the  experimental  data  below  10  keV  represents  a  best 

guess  since  this  data  was  only  taken  accurately  on  spherical  targets.^ 
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For  the  1.5  x  10  W/cm  case,  effects  of  lateral  conduction  are  shown  for 


3  different  conduction  lengths  Lc  -  15p,  50y  and  infinite  conduction 
length  (no  loss  is  assumed  in  the  radial  direction).  For  the  other  cases, 
no  radial  loss  is  modeled.  For  the  5  x  10  W/cm  case,  two  numerical 
x-ray  spectra  are  shown.  That  one  without  a  hard  x-ray  tail  has  been 
obtained  by  turning  off  ion-acoustic  turbulence.  The  lack  of  a  hard  x-ray 
tail  in  that  case  can  be  directly  related  to  the  reduced  heat  flux  inhi¬ 
bition  since  the  extra  absorption  due  to  ion  acoustic  turbulence  is  quite 
small  (39%  total  absorption  compared  to  32%  when  the  ion  acoustic 
turbulence  is  turned  off). 

As  apparent  in  the  various  figures,  the  agreement  is  quite  satisfactory 
since  the  difference  between  numerical  and  experimental  results  are  at  most 
a  factor  of  10  over  6  orders  of  magnitude.  Moreover,  only  a  small  fraction 
of  the  absorbed  laser  energy  can  be  found  in  the  X  rays.  Taking  into 
consideration  the  fact  that  only  the  input  energy,  pulse  shape  and  target 
material  are  specified  in  the  numerical  model,  without  any  assumption  on 
how  much  energy  is  absorbed  and/or  reflected,  we  find  the  agreement  quite 
satisfactory.  The  hot  electron  temperature  measured  from  the  x-ray  spectra 
between  20  and  90  keV  is  shown  in  Fig.  5  for  the  experiments  and  for  our 
calculations. 

The  remarks  which  can  be  made  on  the  differences  between  the  theoretical 
and  experimental  curves  are: 

-  The  numerical  curves  are  smoother  than  the  experimental  curves.  We 
have  already  brought  up  this  point  in  the  introduction. 

-  The  numerical  x-ray  curves  lie  above  the  experimental  ones  for  the 
low- ir rad lance  case  whereas  they  lie  below  for  the  higher  irradiance  case. 

It  could  be  that  suprathermals  near  the  critical  density  which  increase  in 


number  with  irradiance  might  be  neceasary  to  explain  the  high-irradiance 
case.  Another  possible  explanation  is  that  the  ion  acoustic  fluctuation 
level  which  has  been  assumed  constant  does  in  fact  increase  with  laser 
irradiance,  reducing  the  electron  thermal  heat  flux  and  raising  the 
electron  temperature. 

IV.  Conclusion 

In  this  work,  we  have  shown  that  high-energy  x-ray  spectra  can  be 
explained  with  a  one-temperature  fluid  model  which  incorporates  anomalous 
transport.  The  agreement  is  quite  satisfactory  over  6  orders  of  magnitude 
in  x-ray  intensity  and  over  a  factor  of  6  in  laser  irradiance. 

An  essential  ingredient  of  the  fluid  theory  is  the  presence  of 
anomalous  collisions  which  reduce  the  mean  free  path  of  the  electrons 
from  their  free-streaming  value  down  to  a  few  microns  which  are  charac¬ 
teristic  of  the  scale  length  in  Nd- laser  produced  plasmas. 
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X-RAY  ENERGY  (keV) 

Fig.  1  -  Comparison  between  experimental  (solid  line)  and  numerical 
(dotted  line)  hard  X-ray  spectra  at  1.5  x  10^  W/cm  .  The  various 
numerical  spectra  correspond  to  different  value  of  the  lateral  con¬ 
duction  length  L  . 
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Fig.  2  -  Comparison  between  experimental  (solid  line)  and  numerical 
(dotted  line)  hard  X-ray  spectra  at  2.5  x  10^  W/cni  .  The  numerical 
spectrum  has  been  obtained  with  no  radial  loss.  Note  that  for  this 
case,  the  numerical  curve  lies  above  the  experimental  one. 
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Fig.  3  -  Comparison  between  hard  X-ray  spectra  at  5  x  10  w/cm  .  Note 
the  very  good  agreement  between  the  experimental  (solid  line)  and  the 
numerical  (dotted  line)  spectra.  The  numerical  curve  with  the  ion 
acoustic  turbulence  turned  off  points  out  the  importance  of  heat-flux 
limitation. 
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Fig.  4  -  Comparison  between  h 
numerical  curve  (dotted  line) 
line) . 


